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Abstract 

We study the collider phenomenology of the leptophobic Z' boson from an extra U{1)' gauge symmetry in 
models with N-Higgs doublet fields. We assume that the Z' boson at tree level has (i) no Z-Z' mixing, (ii) 
no interaction with the charged leptons, and (iii) no flavour-changing neutral current. Under such a setup, 
it is shown that in the iV = 1 case, all the U{1)' charges of left-handed quark doublets and right-handed up- 
and down-type quarks are required to be the same, while in the N > Z case one can take different charges 
for the three types of quarks. The N = 2 case is not well-defined under the above three requirements. We 
study the pp — Z'V —>■ hhV processes (U = 7, Z and lU^) with the leptonic decays of Z and at the 
LHC. The most promising discovery channel or the most stringent constraint on the U(l)' gauge coupling 
constant comes from the Z'j process below the tt threshold and from the tt process above the threshold. 
Assuming the collision energy of 8 TeV and integrated luminosity of 19.6 fb“^, we find that the constraint 
from the Z''y search in the lower mass regime can be stronger than that from the UA2 experiment. In the 
iV > 3 case, we consider four benchmark points for the Z' couplings with quarks. If such a Z' is discovered, 
a careful comparison between the Z'7 and Z'W signals is crucial to reveal the nature of Z' couplings with 
quarks. We also present the discovery reach of the Z' boson at the 14-TeV LHC in both iV = 1 and N > 3 
cases. 
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I. INTRODUCTION 


After the LHC Run-I completed its operation, we have learned that there exists a scalar boson 


with the mass of about 125 GeV and properties consistent with the Higgs boson in the standard 
model (SM) [l, This suggests that the Higgs sector consist of at least one isospin doublet 
scalar field. At the same time, there has been no report about any other new particles, imposing 
constraints on the parameter space of new physics models, particularly their masses and couplings 
with the SM particles. 

Among new physics models, models with an additional broken U{iy gauge symmetry provide 
one of the simplest framework and have been discussed based on various motivations, e.g., grand 
unified theory (GUT) models j^. This class of models features in an extra massive neutral gauge 
boson Z' whose properties strongly depend on the nature of the U{1)' symmetry. Therefore, 
detection of the Z' boson and detailed measurements of its properties would be a direct probe of 
new physics beyond the SM. 

Searches for such Z' bosons have been done in various collider experiments. The golden search 
channel for models with significant Z' couplings to charged leptons is the Drell-Yan (DY) process 
with e~^e~/fj,~ final states. For example, if the Z' couplings to the SM fermions are exactly the 
same as those of the SM Z boson, i.e., the so-called sequential Z' scenario, the Z' mass mz' is 
constrained to be larger than 2.86 TeV at 95% GL using the Z' —)• channels at the 

LHG with the collision energy of 8 TeV and integrated luminosity of 19.5 fb“^ j^j. 

However, if Z' does not couple or couples very weakly with the charged leptons, the DY channel 
with leptonic final states is no longer useful and one has to resort to hadronic channels. Such a 
leptophobic Z' boson^ can be realized in some GUT models, as is well-known in the Eq model Q]. 
In this case, the most stringent experimental bound on the mass has been given by the data of 
pp —)• ti process at the LHG when mz' is above the tt threshold, with the exact value of extracted 
lower bound depending on the scenario. For example, the lower mass bound is about 1 TeV in 
Scenario I of the Eq GUT model defined in Ref. [g]. 

In our previous work we have studied the LHG phenomenology of the leptophobic Z' boson 
inspired from the Eq GUT model, concentrating on the mass regime below the ti threshold. There 
all the Z' conplings to quarks are determined uniquely according to a given embedding scheme of 


^ The collider phenomenology of a leptophobic Z' boson which also couples to dark matter has been discussed in 
Ref. 01 . 
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Eq However, the Higgs sector, the Yukawa couplings, and their consistency with the leptophobic 


scenario were not discussed in that work. Clearly, more than one Higgs doublet field is required to 
construct the Yukawa Lagrangian that provides the masses of all the quarks and charged leptons 
if different U{iy charges are to be assigned to fermions [l^ . 

Since a SM-like Higgs boson has been discovered, we would like to examine in this work how 
the leptophobic condition affects the structure of the Yukawa sector and study the corresponding 
collider phenomenology. More specihcally, we consider models having a leptophobic Z' boson 
associated with a 1/(1)' gauge symmetry and N Higgs doublets charged under the new symmetry. 
In order to be phenomenologically viable, we require that the Z' boson at tree level have (i) 
no Z-Z' mixing, (ii) no interaction with the charged leptons, and (iii) no flavour-changing neutral 
current (FCNC). With these conditions, we hrst derive consequences about the couplings of the 
leptophobic Z' boson with quarks. Next, we consider constraints on the Z' mass m^i and the U (1)' 
charges of the quarks by using current data. For the mass regimes of niz' < 2rm and mz' > 2mt, we 
take into account the constraints from pp ^ Z' ^ jj at the UA2 experiment [l^ and pp ^ Z' ^ tt 
at the LHC 17|], respectively. Finally, we perform simulations of the pp —)• Z'V —)■ bbV (Y = 7 , Z 
and W^) processes with Z and decaying leptonically at the LHC for N = 1 and N > 3 cases. 

The structure of this paper is as follows. In the next section, we write down the Lagrangian 
satisfying the above-mentioned three conditions, focusing on the Yukawa sector and deriving the 
Z' couplings with the quarks. The Z' —)■ bh decay branching ratio is also derived for later uses. 
In Sec. HI, we discuss the constraints on the 11(1)' gauge coupling constant as a function of the 
Z' mass from the UA2 experiment and the pp ^ Z' ^ tt process at the LHC. Sec. IV shows 
detailed simulations of the t-channel pp —)• Z'V —)• bbV processes at the LHC. We explicitly study 
the one-Higgs doublet scenario and four benchmark points for the N > 3 case. We summarize our 
hndings in Sec. V. 


^ There are six phenomenologically distinct schemes in the Eq model [l5 . [lll| , all of which are examined in Ref. . 
® In general, the U{1)' gauge anomaly cannot be canceled within the SM particle content. Presumably, the anomaly 
is canceled by introducing new particles with non-zero U{lY charges at higher mass scales, which are not discussed 
in this work. 
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II. THE LEPTOPHOBIC Z' BOSON 
A. Lagrangian 


Consider the model with a leptophobic Z' boson associated with a broken U{1)' gauge symmetry 
and N Higgs doublet helds, all assumed to participate in electroweak symmetry breaking. The 
most general kinetic terms, including kinetic mixing with the angle interaction terms for Higgs 
doublet fields (i = 1,... ,N) and interaction terms for a fermion ijj are given by, respectively, 




1 sin X 
siny 1 


N 


= Y^li-igT-W- - ig'YB^ - i~gz'Q' 


2=1 


c-tt = + g'YB^ + • 



with X^i, 


d^x,-d,x^, ( 1 ) 

( 2 ) 

(3) 


In the above Lagrangian, kL“, B^j^ and Z'^ are the gauge fields for the SU{2)l, U{1)y and t/(l)' 
gauge groups, respectively, and the corresponding coupling constants (generators) are denoted by 
g (T“), g' (y) and gz' {Q')- Through a non-unitary transformation. 



(4) 


the mixing term in vanishes in the basis of Z'^ and However, they can still mix with each 
other through terms in after the Higgs fields develop vacuum expectation values (VEV’s), as 
discussed below. 

In the new basis, the interaction terms are rewritten as 

N 


’b _ 


where gz' = gz’ / cos x, and 


Y^\{-igT^W; - ig'YB^ - igz^QZ'^)^^ , 

2=1 

(5) 

+ g'YB^ + gz'QZ'^)^ , 

(6) 

\= Q' + 6Y , with 6 = — tan x ■ 

gz’ 

(7) 
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From Eq. ([5|), the mass terms for the neutral gauge bosons are 





-gg'v^ 

-2ggz' Y.i v‘jQ{^i)^ 



1 

to 

g''^v^ 

2g'gz' Tji 


B^^ 

\-2ggz' Ef 2g'gz> Ef 

44'Ef ) 




( 






).A / z^\ 

\Z'^) 


( 8 ) 


0 -2gz9z' Ef vfQ{<^>i 
0 0 0 
\-29z9z'Ef'vfQm 0 ) 

where Vi {i = 1,... N) are the VEV’s of <f>j and satisfy the sum rule Ei u? = = {V2Gp)~^, with 

Gp being the Eermi decay constant and = sjg^ A g'^ ■ In the last step of the above expression, 
the SU{2)p and U{1)y gauge fields are rotated to the mass eigenbasis in the usual way: 


''V 
B,. 


cos Ow sin 6w 
— sin 6\y cos 6w 



(9) 


where 9w is the weak mixing angle. 

In this paper, we restrict our considerations to the Z' boson that at tree level has 


(I) no mixing with the Z boson; 


(E) no interactions with left- and right-handed leptons (leptophobic condition); and 
(in) no ECNC via neutral scalar bosons 

In Eq. ([8]), the term Ei^i gives rise to non-zero mixing between Z and Z'. We therefore 

impose 

N 

Y.vfQ{^i) = o ( 10 ) 

i=l 

to satisfy (I)^. Secondly, the leptophobic condition (E) demands 

Q{Ll) = Q{eR) = 0 , (11) 

where Lp and ep are respectively the left-handed lepton doublet field and the right-handed charged 
lepton field. Einally, we consider condition (El). In general, a model with a multi-doublet struc¬ 
ture has ECNC’s at tree level because the fermion mass matrix may not be proportional to the 


Tree-level FCNC’s via neutral gauge bosons are automatically forbidden by the Glashow-Iliopoulos-Maiani mech¬ 
anism [28l| . 

® This requirement can be guaranteed when the Z' boson does not couple with the Higgs boson, but obtains its 
mass from the VEV of another scalar boson that breaks the 11(1)' gauge symmetry. 


5 












corresponding Yukawa interaction matrix. In that case, the interaction matrix is non-diagonal in 
the fermion mass eigenbasis. To avoid such a situation, we require that each of up-type quarks, 
down-type quarks, and charged leptons couple to only one Higgs doublet; namely, the Yukawa 
Lagrangian in the N Higgs doublet model assumes the following form; 

= —YdQL^ddR — YuQl^uUR — YeLi^eBR + h.C., (12) 

where and are same or different Higgs fields oi Qr, dj^ and Uj^ represent respectively 

the left-handed quark doublet field, the right-handed down-type quark field and the right-handed 
up-type quark field. For the Lagrangian in Eq. (fT^ to satisfy condition (M), should be all 

distinct; i.e., 

for , (13) 

because both <I>j and could couple to the same type of fermions otherwise. According to Eq. ()12p . 
the Z' charges of the particles have the relations: 

- Q'{Ll) + Q'(ch,) + Q'icR) = -Q{Ll) + Q($e) + Q(eij) = 0 , (14) 

- Q!{Ql) + Q\^d) + Q'{dR) = -Q{Ql) + Q($d) + Q{dR) = 0 , (15) 

- Q!{Ql) - Q'i^u) + Q'iuR) = -Q{Ql) - ^(‘h^) + Q{ur) = 0 , (16) 

where Eq. ([7|) has been used. From Eqs. (fTT]l and (fTTl) . we have 

Q{^e) = 0 . (17) 

That is, the Higgs doublet that couples to the charged leptons cannot carry nonzero Z' charge. 

We now discuss consequences of Eqs. (fT5]) and (fT 6 ]l for several special cases of N. li N = 1 , 
corresponding to the case where <I>e = 4*^ = $„, the Z' charges of all quarks have to be the same: 

Q{Ql) = Q{ur) = Q{dR) . (18) 

In the case of 77 = 2 , the conditions in Eqs. (fTOll and (fT71) require the Z' charges of the two doublets 
to be zero. However, this is in contradiction with Eq. (na). Therefore, the two-doublet case cannot 
simultaneously satisfy all the requirements (I), (E) and (HI). 

In the case of > 3, we obtain from Eqs. (|15p and (|16l) that 


Qi^d) = Q{Ql) - Q{dR) , 

(19) 

Q(^>n) = -QiQl) + Q{ur) . 

(20) 
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In addition, from Eqs. m and (I17p we have the relationship among the Higgs VEV’s and Q 
charges as 


i^u,d 


( 21 ) 


where the sum is over all Higgs doublets other than and <I>rf that do not couple to fermions. In 
the case of = 3 in particular, Eq. (|2ip can be rewritten as 

2 


Q(<I>rf) /(chO)Y 
Qi^u) 


( 22 ) 


Therefore, the three-Higgs doublet case is the minimal setup that allows different Z' charges for 
all the three types of quark fields; i.e., Q{Ql), Q{uji), and Q{dj^) being all different. Moreover, 
the two charges Q(<hrf) and Q(<I>u) have opposite signs according to Eq. ([2^ . Also, it should be 
noted that if (‘h°) and ($[)) are the same, Q{uf^) = Q{d^ according to Eqs. (fT9l) and (M. 

Even though we have explicitly imposed the condition of Eq. cnD, Z'-Z mixing [sj, l9| can still 
occur at loop levels. In our scenario, there are SM quark loop contributions to the two point 
function of Z-Z' mixing at the one-loop level. Its transverse part is calculated as 


nzz'(P^) = 




-E 


r 

/■I 

J „SM 2 

< az^qCLqmq 

/ dx In Ajr — A 

Jo 


+ P^{vz^Vq + af^^qag) 


SM, 


f (ix(x^ — 2x + 1/2) In Ai?— 

Jo 


A' 

"6 


(23) 


where (®|^) Vq (aq) are respectively the vector (axial-vector) couplings of the Zqq and 
Z'qq vertices with corresponding expressions given in Eqs. (1271) and (1721) . Ap = -x(l - x)p^ - niq, 
and A is the divergent part of the loop integral. In the MS scheme, A is simply replaced by In 
with /i being an arbitrary mass scale. Such one-loop contributions give rise to nonzero off-diagonal 
(1,3) and (3,1) elements of the mass matrix Eq. ([8]). In this case, the Z-Z' mixing can be calculated 
for given values of momentum and scale /r as 


tan29zz' — 


2nzzJp‘^) , [nzzjp^)? 


2 2 
— m^, 


+ o 




rur 


rur. 


(24) 


Taking = 150 GeV, g^: = 0.1, Vu = Vd = 1 {vu = Vd = 0.5) and a„ = = 0 (a^ = ad = 0.5), 

we obtain sin Ozz' = 1-65 (1.04) x 10“^ with = gj = rin? 2 :i- this calculation, we use mt = 173.07 
GeV and mf, = 3.0 GeV, and all the other quark masses are neglected. Although such a mixing 
effect can contribute to additional Z' boson productions in collider experiments, it is negligibly 
small because 9zz' = 0(10“^). 
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B. Z' couplings to quarks 


In a model with N Higgs doublet fields, the Z' interactions with one generation of quarks are 
given by Eq. Q as 


^fnt = 9z' Q{Ql)Qi^PlQ + Q{uR)u'y^PRU + Q{dR)d'y^PRd 




(25) 


where the projection operators Pl,r = (1 T 75)/2- Alternatively, Eq. (1251) can be written in terms 
of the vector coupling Vq and axial-vector coupling Uq as 


Pint = 9z'qi^{vq - l5aq)qZ'^ , 


(26) 


where 

^ [Q(Ql) + Q{qR)] , flg = ^ [2 (Ql) - Qign)] , for g = u, d. (27) 

When Ql carries a nonzero Q charge, one can always normalize its value to unity by rescaling the 
coupling gz'- Therefore, we take the value of Q{Ql) to be either 1 or 0 in the following analyses. 
In this convention, aq is equal to 1 — {—Vq) when Q{Ql) = 1 (0). In the one-Higgs doublet case, 

only gz' is a free parameter and all the others are fixed according to Eq. m as: 


Vu = Vd = l, au = ad = 0 . (28) 

This means that the Z' couplings to the quarks must be vectorial. The interaction Lagrangian in 
Eq. (1261) can also be rewritten in terms of the chiral couplings as 

Pit = 9z'qi>^{9l Pl + 4 PR)qZ'^ , (29) 

where 


gl = Vq + aq = 1 , g‘}^ = Vq-aq = 2vq-l . 


(30) 


C. Z' Decays 


The partial width of the Z' decaying into a pair of quarks is given by 


T{Z' qq) = 9z>^^ [^^^(1 + 2xg) + 0^(1 - 4xg)] - Axq , with Xq = 


m„ 


m 


z' 


The corresponding branching fraction is 


Br(Z' qq) = 


r{Z' ^ qq) 
Ez' 


(31) 


(32) 








Br{T^bb) 



-0.5 o!o' ' 0^5 ' ' LO L5 

Vu 


FIG. 1: Contour plot of the Z' —>• hh branching fraction on the Vu-Vd plane, assuming mz' < 2mt and 
ignoring the quark masses. 


where the total decay width Tz> = QQ) summed over all quarks with mass less than 

mz 'When mz' < 2mf but much greater than 27T7-t,, the branching fraction is approximately 


Br(Z' —^ qq) 


2{vg - 1 / 2 )^ + 1/2 

6 (urf - 1/2)2+ 4(u„-1/2)2+ 5/2 ’ 


(33) 


where Uq = 1 — Vg is used and the quark masses have been neglected. Among the various decay 
modes, the Z' —)• bb decay channel with b-tagging can be the most important one for discovering the 
Z' boson at colliders, especially in the small mass regime mz' < 2mt. Fig. [T]shows the contour plot 
of Br(Z' —)• bb) on the v^-Vd plane. The branching ratio increases (decreases) as Vd {vu) deviates 
from 1/2. 


III. CONSTRAINTS 


We first consider the constraint on the coupling constants of the Z' boson to quarks from the 
UA2 experiment. The UA2 experiment had searched for a Z' boson via the pp —)• Z'^*'^ —)• jj 
process at the center-of-ma ss ( CM) energy of 630 GeV. The analysis was done in the mass range 
from 100 GeV to 300 GeV [l3|. The cross section for the hard process qq —)• Z'^*^ —)■ q'q' is given 
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FIG. 2: TZpp as a function of raz' at the CM energy of 630 GeV using the CTEQ6L PDF’s. 


by 


^qq^q'q' 


{s) = 


9z' 


(vq + a^n){vi + ai) 


(34) 


127rs (1 — m?^i/sY + ’ 

where \/I is the CM energy of the partons q and q. As can be inferred from Eq. (|3ip . Tz' I’mz' ~ 
g'|,/(47r) whose value for our benchmark points defined below is of 0(0.01). Therefore, the narrow 
width approximation can be employed to simplify the cross section as 


^qq—>-q'q‘ 


q'{s) ^ 


^- ^(Vq + ag)iVq' + a^')^(s - ml,) 

12 m^,l z' ^ ^ ^ 

+ Qq) X Br(Z' q'q')6{s - m|,). 


(35) 


The cross section for the pp —)• Z'* —> q'q' process is obtained by convoluting the above expression 
with the partonic luminosity functions dCqq/dT{T, pp) for the qq initial state as 


/■l pr 

a{pp q'q') ='^ dT-^{T, PF)dqq^q'q'{s = Ts), 




(36) 


where fip is the factorization scale, and ^/s = 630 GeV. The luminosity function is given by 

d^ggi'^,9F)= - fq{x,pF)fq{T/x,pF), (37) 

Jt ^ 


with fq and fq being the parton distribution functions (PDF’s) of q and q, respectively. Again, in 
the narrow width approximation, the cross section becomes 


a{pp Z' ^ q'q') ~ '^Cqq{ml, / S, PF)dqq^q'q'{s = m|,) 


5 •-dd 


g,g 


( Vu 


1 \^ / 1 \^ 1 


+ \Vd - -j + -^{T^pp + 1 ) 


Br(Z' ^ q'q'), 


(38) 
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FIG. 3: Contours of the pp ^ Z' ^ jj cross section with the CM energy of 630 GeV on the Vu-Vd plane, 
where gz' = 0.18 (0.26) and mz' = 150 (300) GeV are applied to the left (right) plot. The thick black 
curves correspond to the maximum cross section allowed by the UA2 experiment. 

where TZpp = Cuul^^dd denotes the ratio of luminosities of the uu and dd initial states. It is seen 
that the cross section is proportional to the elliptical expression of Vu and Vd that has dependence 
on TZpp. In Fig. [2l we show TZpp as a function of mz' in the mass of interest to us. Here fip = w-zo 
^/s = 630, and the CTEQ6L PDF’s are used. 

According to Eq. (I38p . the maximum of (vd) for given values of gz' and the cross section 
is obtained when we hx vp (vu) = 1/2. In particular, the same cross section for {vu,Vd)={l,l) is 
obtained by taking {vu,Vd)={l/2,v'^^^) and {vu,Vd)={v^^^,1/2) for the same value of where 

<■“ = + I’ (fo-- = 1/2), (39) 

vT"" = + + (forn„ = l/2). (40) 

In Fig. [3l we show the contours of the pp ^ Z' ^ jj cross section for mz' = 150 GeV and 
gz' = 0.18 {mz' = 300 GeV and gz' = 0.26) in the left (right) plot. These values of gz' are 
obtained for the choice of (?;„,'(;rf)=(l,l). The thick contour in the left (right) plot corresponds to 
a ~ 150 (5.2) pb, the upper limit set by the UA2 experiment {l^ . From the thick contours, we 
obtain the values of to be about 1.8 and 2.5 in the left and right plots, respectively, as can 
also be obtained by using Eq. (jiOl) with TZpp = 5.4 and 14.3 given in Eig. [2l 

In the following, we estimate the upper limit on the gauge coupling g^, using the UA2 data 
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Vu 

Vd 

„max 

yz' 



BPl 

1 

1 

max(BPl) 

9z' 

0 

0 

BP2 

1/2 

„,max 

max(BP2) 

9z' 

1/2 

1 - 

BPS 

1/2 

1/2 

max(BP3) 

9z' 

1/2 

1/2 

BP4 

0 

1/2 

max(BP4) 

9z' 

0 

-1/2 


TABLE I: Four benchmark points defined by the values of their vector couplings, and Vd, and the 
maximum allowed by the UA2 experiment. BPl, BP2 and BPS have Q{Ql) = 1, while BP4 has 
Q{Ql) = 0. The value of is a function of mz' shown in Fig. |4l Also shown are the corresponding 
axial-vector couplings a„ and ad for quick reference. 



FIG. 4: The upper limit of the coupling of Z' for several benchmark points which are obtained from the 
Pp jj cross section measurement of the UA2 experiment. 

for several benchmark points (BP’s) of {vu,Vd) defined in Table HI BPl corresponds to the one- 
Higgs doublet case. BP2 has the maximum value for Vd using , the upper limit of g^i for 

BPl. BPS corresponds to the case with purely left-handed Z' couplings; i.e., given in Eq. (ISOl) 
vanishes. BP4 is the case with nonzero Z' couplings only for the right-handed down-type quarks; 
i-e., Q{Ql) = Q{ur) = 0. 

We compute the cross section of the pp ^ Z' ^ jj process at the CM energy of 630 GeV using 
MADGRAPH/MADEVENT 5 [l^ and our model files, and find the maximum coupling for each 

benchmark point that saturates the cross section upper bound at 90% confidence level (CL) shown 
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FIG. 5: Feynman diagrams for the gauge boson associated production process of the Z' boson at the LHC. 


in Fig. 2 of Ref. [13|]- In this calculation, we take the Z' width calculated by CalcHEP3.6.15 [R 
We also apply a global K-factor K = 1.3 for the cross section [l^. Fig. [J] shows the upper limit of 
the gauge coupling for each BP, where the limit for BP2 is taken to be the same as that for BPl, 
and that for BP4 is consistent with the constraint given in Fig. 1 of Ref. [l^. We will use these 
BP’s and the corresponding constraints in the following studies of collider phenomenology. 

In the case of mz' > 2mt, the s-channel pp —)■ Z' —)• ti process is most useful for the Z' search 
at the LHC. The CMS group reported the search for production of heavy resonances decaying into 
tt pairs using the data of integrated luminosity of 19.6 fb“^ at 8 TeV IT]]. Nonobservation of an 


excess in this process provided an upper limit on the cross section (pp —)• Z') times the branching 
fraction of Z' —)• tt at 95% CL as a function of the invariant mass M^i of the tt pair. Comparing this 
limit with the cross sections of pp ^ Z' ^ tt computed for our scenarios, one can also obtain the 
constraint on gz' as a function of mz'- In our calculation, we apply a global K-factor of iL = 1.4 
taken from Ref. 1^. This bound in the mz' > 2mt regime will be imposed on the study of each 


BP in the next section. 


IV. COLLIDER PHENOMENOLOGY 


A. Gauge boson associated production of Z' 


Consider the t-channel production of the Z' boson associated with a gauge boson V (= 7 , Z or 
W^) at the LHC, i.e., the pp —)• Z'V processes. The dominant contributions of the Z'V production 
are given by quark initial states qq {qq') —)• Z'V as shown in Fig. O For the Z' 7 , Z'Z final states, 
additional contributions come from the gluon fusion processes through quark box diagrams. Their 
cross sections are only a few percent of the dominant processes at the LHC [^]. We thus neglect 
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FIG. 6: Contour plots for the cross section of the pp —>■ Z'j (left) and pp —>■ Z'Z (right) processes for 
the 8 -TeV LHC on the Vu-Vd plane, where gz' = 0.18 and niz' = 150 GeV are applied. The black curves 
correspond to those in Fig. |3l 

the gluon fusion contributions in the following analysis. The pp —)• Z'V cross section is given by 
t7(pp Z'V) (X ^ ^ ^ 0 

+ + 4!i)^ + (41) 

where the SM vector and axial-vector couplings 

4 “) =(Q., 0 ), for y = 7 , 

(4!^, 4!^) = (4/2 - Qq sin^ Ow, 4/2) , for F = Z, 

(4!^i4!^) = for y = iy, (42) 

with Qq and Iq being respectively the electric charge and the third isospin component of the quark 
q. In Eq. (HH), TZpp is the ratio of the luminosity functions for the pp collision, analogous to TZpp 
in Eq. (1381) . Eor the collision energy of 8 TeV and m^, = 150 (300) GeV, TZpp is found to be about 
1.5 (1.6) using the CTEQ 6 L PDE’s. 

Fig. [ 6 ] shows the contour plots of the cross section for the Z'’j (left) and Z'Z (right) final states 
on the {vu,Vd) plane, where we take gz' = 0.18 and mz' = 150 GeV as in Fig. [3l Gontours of the 
pp ^ Z' ^ jj cross sections in Fig. [3] are also shown in the plots for comparison. The pp —)• Z'W 
cross section is constant, about 5.5 pb in the case of gz' = 0.18 and mz’ = 150 GeV, on the {vu,Vd) 
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plane as shown by Eqs. m and (H2|) . It is observed that the cross sections of the pp —)• Z'-y and 
pp —)• Z'Z processes have similar dependence on {vu,Vd) as that of the pp ^ Z' ^ jj process. This 
can be readily understood as follows: the dependence on {vu-,Vd) is determined by the elliptical 
expression given in Eq. (|41l) similar to the cross section of pp —>• Z' —)• jj. The shape of ellipse is 
determined by [Qu/Qd^T^pp according to Eqs. ()4T]1 and (H2D . Eor the case of m^, = 150 GeV, this 
factor is about 6 , close to TZpp. 

As mentioned in Section III Cl the Z' —)• bb decay is the most promising mode to search for Z' 
in the low-mass regime as one can use b-tagging to reduce backgrounds. We consider the signal 
events 


pp 

pp 

pp 


Z'7 667 jbjb'y, 

Z'Z bbi+r jbjb^^r, 
Z'W^ ^ bbi^^T^ jbjbi^^T, 


(43) 

(44) 

(45) 


where i = e ov p, .^-ris the missing transverse energy, and jb is a tagged b-jet. For these production 
processes, we also take into account the enhancement from QCD corrections characterised by the 
K-factor K = 1.3 j^, 23| in our analysis. The SM backgrounds corresponding to each of above 
signals are 


pp 667 jbjbi, 

PP -t JJT jbjbl, 

(46) 

pp bbi+£~ jbjb^^^~, 

PP -t jj£^£~ -t jbjbi^^~, 

(47) 

pp bbi^^T -t jbjb£^^T, PP -t jj£^^T -t jbjb£^^T, 

(48) 


where j is a jet coming from a gluon or a non-bottom quark. Since the backgrounds involve various 
processes with different K-factors, some of which have not been evaluated yet, we take AT = 1.2 
and 1.4 to estimate possible uncertainties. The signal and background events are both generated 
using MADGRAPH/MADEVENT5, and passed to PYTHIA6 via the PYTHIA-PGS package to include 
initial-state radiation, final-state radiation and hadronization effects. We note in passing that in 
MADGRAPH/MADEVENT 5 the factorization and the renormalization scales are set as V!?_;^(M?+p^„)/2 


for the two hnal-state particles. The detector-level simulation is carried out using PGS4 


25l |. which 


performs the b-tagging with an efficiency about 0.5 for high-energy jets in the central region dehned 
by the jet rapidity limit \'r]{j)\ < 2.0. The number of signal events is reduced to ~ 10% due to 
double b-tagging and the rapidity cut. 

In addition, since the b-jets are energetic and boosted along the direction of Z', we thus impose 
the following kinematic cuts for the transverse momentum of each b-jet, PTijb)^ and the rapidity 
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difference between the two b-jets, 


PrUb) > 40 GeV, \Ar]j^jJ < 2.0, (49) 

where the lower and upper limit on PT{jb) and are chosen by optimizing the cut efficiency 

at mz’ = 150 GeV. For the jhjbl events, we further eliminate soft photons by applying the following 
Pt cut: 


Pt{i) > 10 GeV for jbjby , 


(50) 


where we have not chosen a larger lower limit for PTipf) as the photon energy in Z'^ production 
process tends to be small. We also take the following cuts for jbjb^^^~ and jbjb^^^T events: 


PT{i) > 25 GeV for jbjb^~^^ and jbjb^'^^T , 

> 25 GeV for jbjb^-'^^T , (51) 


where the lower limit for pt{^) is taken from the W^bb search at GMS [2^ and the same value is 
used for Finally, we also make a cut on the invariant mass of the two b-jets 


mz'{l — 0.2) < < mz' + 10 GeV . (52) 

Here we have chosen asymmetric limits, as also used in our previous paper Q], because the shape 
of bb invariant mass distribution is not symmetric around mz'- 


B. One-Higgs doublet case 


We hrst apply the above analysis to the one-Higgs doublet case (or BPl), where = 1. 

Here we remind the reader that each BP by definition takes the maximum gauge coupling, 
that saturates the UA2 bound. Therefore, the cross sections given in the following numerical 
calculations are their maxima derived from the curves given in Fig. HI In Fig. [71 the cross sections 
of the pp —>• Z'V processes at the 8-TeV (left plot) and 14-TeV (right plot) LHG are given as 
functions of mz', all computed with CalcHEP. It is seen that the pp —)• Z'^ process gives a ~ 4 and 
10 times larger cross section than the pp —)• Z'W and pp —)• Z'Z processes, respectively. 

In Table IIIl we show the cross sections for the signals and the backgrounds given in Eqs. (|43l) - 
(I45l) and in Eqs. (|46]) - (|48]) . respectively, at the collision energy of 8 TeV. We take mz’ = 150 GeV 
as an example and apply the corresponding upper limit of gz' in Fig. HI The signal significance 


defined as 


27| 


S = v^2[(s -P b) ln(s/6) - s]. 


(53) 
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FIG. 7: Cross sections of the pp — Z'V {V = and Z) processes for the one-Higgs doublet case (BPl) 
at the 8-TeV (left) and 14-TeV (right) LHC. 


Events 

jbjbj 

jbjbi+i- 

jbjbi^Ex 


as [pb] aB° [pb] 5 

as [pb] aB° [pb] S 

as [pb] aB° [pb] S 

b-tagging 

1.6x10"^ 18. 4.9 (4.6) 

1.1x10"® 1.1x10-1 0.43 (0.40) 

1.3x10-2 5.9x10-1 2.2 (2.1) 

Prijb) > 40 GeV 

1.0x10"^ 4.8 5.9 (5.4) 

6.5 xlO""* 2.8x10"^ 0.49 (0.45) 

7.7x10"® 1.9x10-1 2.2 (2.1) 

l^Vjbib 1 2.0 

1.0x10-1 4.5 6.0 (5.5) 

6.3x10-1 2.6x10-2 0.50 (0.47) 

7.6x10"® 1.8x10-1 2.3 (2.1) 

PT{e) > 25 GeV 


4.4x10-1 1.1x10-2 0.52 (0.48) 

6.3x10"® l.lxlO"! 2.4 (2.2) 

> 25 GeV 



4.9x10"® 6.3x10-2 2.5 (2.3) 

Mjbib cut 

6.7x10"^ 1.6 6.9 (6.4) 

3.0x10-1 3.4x10"® 0.64 (0.59) 

3.4x10"® 1.8x10-2 3.1 (2.9) 


TABLE II: Cross sections of signals including the K-factor AT = 1.3 [ps) and background processes at leading 
order (cr^®) after sequentially imposing each of the selection cuts shown in the first column in the case of 
mz' = 150 GeV and ^ = 8 TeV. For the significance 5, we take the integrated luminosity of 19.6 fb“^ 
and apply the K-factor for the background cross sections. Values without (within) parentheses correspond 
to the background K-factor of K = 1.2 (1.4). 


is also given in the last column of each final state with the assumption of an integrated luminosity 
of 19.6 fb“^, where s and b denote the numbers of signal and background events, respectively. The 
number without (within) parentheses corresponds to the backgrounds using the K-factor K = 1.2 
(1.4). From the third to last rows, we show the results after sequentially imposing the kinematic 
cuts in the first column, as given in Eqs. ([Ml), ([50l) . ([5T]) . and Eq. ([52]) . The Z''y process has the 
largest significance S due to its largest signal cross section among all. Although the Z'Z process 
has the smallest background cross section, the signal cross section is also highly suppressed due to 
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FIG. 8 : Significances for Z' 7 , Z'W and Z'Z signals at the 8 -TeV LHC with the integrated luminosity of 
19.6 fb“^ (left) and the 14-TeV LHC with the integrated luminosity of 100 fb“^ (right). All kinematic cuts 
in Eqs. (l4^ . (l50l) . (ISTI) and Eq. (l52l) have been imposed. The K-factor of 1.4 for the backgrounds is used. 

the leptonic branching fraction of Z. 

In Fig. [HJ the signal signihcances for the Z'Z and Z'W^ processes are shown as functions 

of mz', assuming the collision energy and the integrated luminosity respectively to be 8 TeV and 
19.6 fb“^ for the left plot and 14 TeV and 100 fb“^ for the right plot. Here all the kinematic cuts 
in Eqs. (I49p - ()52p have been imposed and K = 1.4 is applied to the backgrounds for a conservative 
estimate. The case ol K = 1.2 for the backgrounds would have slightly better significances. Since 
the UA2 upper limit of gz' is used, these significances are the largest values that one can expect. 
For the Z '7 process, 5 > 2 (5 > 5) in the entire mass region of the plot (m^/ < 200 GeV and 
mz' > 260 GeV) at the 8 -TeV LHG, while 5 > 5 is achieved at the 14-TeV LHG. The other 
two processes Z'Z and Z'W^ give smaller values of S as compared to that in the Z '7 process. 
Especially for the Z'Z process, S is smaller than 2 in the entire mass range considered in these 
figures even in the case of 14 TeV and 100 fb“^. 

As shown partly in the left plot of Fig. [9] below, the ti data from the 8 -TeV LHG has a more 
stringent constraint for niz' > 500 GeV. This mode will still be the most promising search channel 
or impose the most stringent constraint for the same mass regime. Between 2mt and 500 GeV, the 
tt constraint is expected to be worse and may be comparable to the Z '7 mode proposed in this 
work. Therefore, we will put our focus on the region of mz' ^ 500 GeV in the following discussions. 
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FIG. 9: The left plot shows the required gauge coupling constant to render the significance of 2, 3, 4 and 5 
at the 8-TeV LHC with the integrated luminosity of 19.6 fb“^ for the Z'y process in the one-Higgs doublet 
case. The right plot shows the corresponding curves for significance of 2 and 5 at the 14-TeV LHC with the 
integrated luminosities of 100 and 300 fb“^. The black curve gives the UA2 and ti constraints on the gauge 
coupling constant. 

The result for Z'^ process given in Fig. [S] can be translated into the required Z' coupling 
to achieve a specific value of S for each value of mz' using the fact that the cross section (number 
of events) is proportional to g'|/. We have applied K = 1.4 to the backgrounds in these estimates. 
The left plot of Fig. [S] shows the contours of the required g^, for 5 = 2, 3, 4 and 5 at the collision 
energy of 8 TeV and integrated luminosity of 19.6 fb“^. The upper limits from the UA2 experiment 
is also shown in this figure. For the region mz' > 500 GeV, the upper bound on g^, is derived from 
the current LHC data of the pp —)• ti process [l^. From the curve for 5 = 2, we obtain a stronger 
constraint on g^, in comparison with that from the UA2 experiment. Similarly, the right plot of 
Fig. 0 shows the required for 5 = 2 and 5 at the collision energy of 14 TeV and integrated 
luminosities of 100 fb“^ and 300 fb“^. The 95% CL upper limit (5 = 2) is well below the limit 
from the UA2 experiment, which is about 0.06 (0.05) at mz' = 150 GeV assuming the integrated 
luminosity of 100 (300) fb“^. The minimal value of g^, for the discovery (5 = 5) is also below the 
UA2 limit, which is about 0.10 (0.08) at mz' = 150 GeV assuming the integrated luminosity of 
100 (300) fb-k 
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mz' [GeV] 

150 

200 

250 

300 


a{Z'-i)xW{Z' bb) [pb] 

4.0 

1.6 

8.1x10-1 

1.7 

BPl 

a{Z'Z)xW{Z' bb) [pb] 

3.9x10-1 

1.8x10-1 

1.0x10-1 

2.4x10-1 


a{Z'W)xBr{Z' ^ bb) [pb] 

1.2 

5.2x10-1 

3.0x10-1 

6.3x10-1 


cr(Z' 7 )xBr(Z' ^ bb) [pb] 

6.1 

2.8 

1.7 

3.9 

BP2 

a{Z' Z)xBt{Z' —>• bb) [pb] 

6.0x10-1 

2.7x10-1 

1.6x10-1 

3.7x10-1 


a{Z'W)xBT{Z' ^ bb) [pb] 

1.8 

7.8x10-1 

4.5x10-1 

9.8x10-1 


a{Z'-f)xBT{Z' bb) [pb] 

4.0 

1.6 

7.0x10-1 

1.7 

BPS 

a{Z'Z)xBi-{Z' ^bb) [pb] 

7.2x10-1 

3.2x10-1 

1.6x10-1 

4.1x10-1 


a{Z'W)xBr{Z' ^ bb) [pb] 

2.4 

1.0 

4.9x10-1 

1.3 


a{Z'"f)xBT{Z' bb) [pb] 

5.7 

2.3 

1.1 

2.0 

BPl 

a{Z'Z)xBi-{Z' ^bb) [pb] 

3.3x10-1 

2.4x10-1 

1.4x10-1 

2.9x10-1 


a{Z'W)xBT{Z' ^ bb) [pb] 

0.0 

0.0 

0.0 

0.0 


TABLE III: Products of the pp —>■ Z'V cross sections and the branching fraction Br(Z' — bb) in units of pb 
for each benchmark points with the collision energy of 8 TeV assuming mz' =150, 200, 250 and 300 GeV. 


C. A^-Higgs doublet case 

For the A^-Higgs doublet case {N > 3), we apply the analysis for the Z' search discussed in 
Section HV Al to the four BP’s listed in Table HI Table Hill shows the products of the pp —)• Z'V cross 
sections and the branching fraction Br(Z' —?■ bh) for several values of mz' assuming the collision 
energy of 8 TeV. Again, the maximum coupling in Fig. 0] is used for each BP. We first discuss 
some general properties of each BP: 

• BPl is exactly same as the one Higgs doublet case discussed in the previous subsection, and 
is shown here for comparison. 

• BP2 has the maximum value for the pp —)■ Z'y —)■ 667 cross section among the four BP’s as 
the branching fraction for Z' —)• bb mode is larger than BPl by ~ 3/2 owing to the larger 
value of Vd- 

• BP3 has the maximum values for the pp —^ Z'W —)• bbW and pp —^ Z'Z —?■ bbZ cross sections 
as a consequence of the purely left-handed couplings. 

• BP4 has nonzero couplings only for the right-handed d-type quarks. Thus, the pp —>• Z'W 
cross section is identically zero. On the other hand, the cross section of pp —)• Z'7 —)• bb'y is 
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close to that in BP2 although the branching fraction for Z' —)• hh is smaller by ~ 2/3. This 
is because the UA2 limit on g^' in BP4 is the largest. 

BP2 is expected to give the largest significance for the pp —)• Z '7 —)• 667 process due to its largest 
cross section. The purely left-handed (right-handed) couplings in BPS (BP4) result in distinctively 
different relative strengths of the production cross sections of Z' 7 , Z'W and Z'Z from the one Higgs 
doublet case (BPl). It is therefore important to compare signal cross sections between three Z'V 
production processes to distinguish the BPs. In the following, we concentrate on BP2 and BPS, 
and estimate the signal significances and the required coupling strength. BP4 gives a significance 
similar to BP2 in the case of Z' 7 , and its significances for other two processes are too small to be 
useful. 

The upper (lower) plots of Fig. HOI show the signal significances of the Z' 7 , Z'Z and Z'W^ pro¬ 
cesses as functions of mz' for BP2 (BPS), assuming the collision energy and integrated luminosity 
of 8 TeV and 19.6 fb“^ (left plots) and 14 TeV and 100 fb“^ (right plots). Here K=1.4 is applied 
to the backgrounds as in the one Higgs-doublet case. For BP2, we obtain larger significances than 
BPl in Fig. [ 8 ] because of the larger branching fraction Br{Z' —)• hh). The hierarchy of significances 
is similar to that in BPl. 5 > 5 for the Z '7 process in the entire mass range for the 8 -TeV case, 
and the other two processes give smaller values of S. For BPS, we obtain a similar significance for 
the Z''^ process as in BPl. However, the significance of the Z'W process is enhanced by the purely 
left-handed couplings to also have similar values as Z '7 for the entire mass range. The plot shows 
that 5 > 2 (5 > 5) in the entire mass range {mz' < 200 GeV and mz' > 250 GeV) in the 8 -TeV 
case for both Z '7 and Z'W processes. The Z'Z process gives a smaller value of S that is less than 
2 (5) in the entire mass range for the 8 -TeV (14-TeV) case. 

Finally, following the analysis for Z'^ in the one-Higgs doublet case, we derive here the required 
Z' gauge coupling 5 ^, to achieve specific values of 5 as a function of raz'- Here we take iF = 1.4 
for the backgrounds as in the one Higgs-doublet case. The upper left (lower left) plot in Figs. [11] 
shows the contours of required g^, for 5 = 2, 3, 4 and 5 at the 8 -TeV LHC with the integrated 
luminosity of 19.6 fb“^ in the scenario of BP2 (BPS). The upper limits from the UA2 experiment 
and the CMS it data are also shown in these plots. For BP2, we find that one can get a stronger 
constraint on as compared to the one-Higgs doublet case in Fig. O For BPS, the constraint 
is weaker but the relative strength between gz' values of UA2 limit and that on fixed-significance 
contours are similar to the one-Higgs doublet case (BPl). This is because BPS produces a ratio 
of the cross sections of pp —>• jj and pp —)• ^Z'{-^ bh) that is similar to the one-Higgs doublet 
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FIG. 10: Same as Fig. |S1 but for BP2 (upper plots) and BPS (lower plots). 

case. The upper right (lower right) plot in Figs. [11] shows the required ( 7 ^, for 5 = 2 and 5 for 
BP2 (BP3), assuming the 14-TeV LHC with the integrated luminosities of 100 fb“^ and 300 fb“^. 
The 95% CL upper limit (5 = 2) is well below the limit from the UA2 experiment, and is about 
0.05 (0.04) and 0.08 (0.07) for BP2 and BP3, respectively, at mz' = 150 GeV and assuming the 
integrated luminosity of 100 (300) fb“^. The minimum value of < 7 ^, for the discovery (5 = 5) is also 
below the UA2 limit, which is about 0.08 (0.06) and 0.14 (0.11) for BP2 and BP3, respectively, at 
niz' = 150 GeV and assuming the integrated luminosity of 100 (300) fb“^. 
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FIG. 11: Same as Fig. HI but for BP2 (upper plots) and BPS (lower plots). 

V. CONCLUSIONS 

We have studied the collider phenomenology of the leptophobic Z' boson from an extra U{1)' 
gauge symmetry in models with multiple Higgs doublet fields under the conditions of (I) no mixing 
with the Z boson, (H) no interactions with charged leptons, and (HI) no FCNC via the Z' mediations. 
We have shown that in the case ol N = 1, all the Q charges for quarks have to be the same; i.e., 
Q{Ql) = Q{ur) = Q{dR)- This consequence is relaxed in models with N > 2. We have explicitly 
shown that the three conditions cannot be simultaneously met in the case of = 2. 

We have discussed the constraint on the 11(1)' gauge coupling constant from currently available 
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data of the UA2 experiment and the CMS pp —)• tt measurements at the LHC. The upper limit 
on the gauge coupling constant is derived for m^/ < 300 GeV from the former and for m^/ > 500 
GeV from the latter. 

We have studied the pp —>• Z'V —)• bbV processes (H = 7 , Z and W^) with the leptonic decays 
of Z and at the LHC for the N = 1 and N > 3 cases. We propose that in the case of = 1 
and mz' less than the tt threshold, the Z'j process serves as the most promising discovery mode 
or provides the most stringent constraint on the U{iy gauge coupling constant, stronger than that 
from the UA2 experiment, at the collision energy and integrated luminosity of 8 TeV and 19.6 
fb“^, respectively. The tt mode is still the best search channel above the tt threshold. 

When N > 3, we have considered four benchmark points (BPl, BP2, BPS and BP4) for the Z' 
couplings with quarks. The benchmark point BPl exactly corresponds to the N = 1 case. In BP2, 
the couplings are chosen such that the Z' —?■ bb has the largest branching ratio. It has a slightly 
stronger constraint on the gauge coupling from the Z'’j process than that in the = 1 case. We 
have also found that the benchmark point with purely left-handed couplings (BPS) gives similar 
values of significances for the Z'j and Z'W processes in contrast to the others. This example shows 
that a careful comparison between the cross sections of Z '7 and Z'W processes will be crucial to 
reveal the nature of Z' couplings with quarks. The scenario of BP4 has nonzero Z' couplings only 
for the right-handed down-type quarks. This results in null cross section for the Z'W process. 
Nevertheless, the Z '7 cross section approaches that in BP2 because of the larger gauge coupling 
strength allowed by UA2. Finally, we have computed the discovery reach of such a Z' boson at 
the 14-TeV LHC in both N = 1 and N > 3 cases. For BPl, the expected 2cj upper limit on the Z' 
coupling g^, has been obtained to be about 0.049 (0.084) for m^, = 130 (600) GeV assuming the 
integrated luminosity of 300 fb“^. In addition, the required Z' coupling to reach the 5 cj discovery 
has been found to be 0.078 (0.13) for = 130 (600) GeV. 
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